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ABSTRACT 

A significant fraction of the hot Jupiters with final circularized orbital periods of < 5 days are 
thought to form through the channel of high-eccentricity migration. Tidal dissipation at successive 
periastron passages removes orbital energy of the planet, which has the potential for changes in semi- 
major axis of a factor of ten to a thousand. In the equilibrium tide approximation we show that, 
in order for high-eccentricity migration to take place, the relative level of tidal dissipation in Jupiter 
analogues must be at least 10 times higher than the upper-limit attributed to the Jupiter-Io interaction. 
While this is not a severe problem for high— e migration, it contradicts the results of several previous 
calculations. We show that these calculations of high— e migration inadvertently over-estimated the 
strength of tidal dissipation by three to four orders of magnitude. These discrepancies were obscured 
by the use of various parameters, such as lag time r, tidal quality factor Q and viscous time ty. We 
provide the values of these parameters required for the Jupiter-Io interaction, tidal circularization and 
high-e migration. Implications for tidal theory as well as models of the inflated radii of hot Jupiters 
are discussed. Though the tidal Q is not, in general, well-defined, we derive a formula for it during 
high-eccentricity migration where Q is approximately constant throughout evolution. 
Subject headings: 



1. INTRODUCTION 

The recent explosion in the discovery of extrasolar 
planets has re-invigorated interest in theory of tidal exci- 
tation and dissipation. The proximity of the hot Jupiters 
to the stellar host suggests that tidal forces shaped or are 
currently shaping the orbital and structural properties of 
these close-in planets. 

High-e migration (e.g., Rasio & Ford 1996; Wu & 
Murray 2003; Fabrycky & Tremaine; see Socrates et al. 
(2011) for a summary) is a class of migration scenar- 
ios that depend upon the excitation of eccentricity, e, to 
high values where tidal dissipation is activated as small 
periastron distances. Successive periastron passages con- 
vert orbital energy into heat that is ultimately radiated 
away. The discovery that the orbits of a significant frac- 
tion of the hot Jupiters are mis-aligned with the spin of 
their stellar hosts lends strength to the claim that some 
high— e migration is at work (Winn et al. 2010). 

For such mechanisms to work, tidal dissipation during 
periastron passage must be strong enough to dissipate 
the required amount of orbital energy within the age of 
the system, which for FGK stars on the main sequence 
is a significant fraction of the Hubble time. This require- 
ment places a lower bound on the strength of dissipation, 
which can then be tested against theories of the tidal ex- 
citation and dissipation. By utilizing equilibrium tidal 
theory in the weak friction approximation, Goldreich & 
Soter (1966; hereafter GS) produced the first estimate of 
the level of dissipation in gas giants (Jupiter and Saturn), 
parameterized by a tidal quality factor Q. 

There is a general understanding in the literature that 
high-e migration is a viable mechanism for forming hot 
Jupiters in the event that the relative level of tidal dissi- 
pation, within the context of equilibrium theory, is con- 
sistent with that inferred by GS's study of the Jupiter-Io 
interaction. Here we critically re-examine the viability 



of high-e migration and tidal circularization for hot ex- 
trasolar gas giant planets. We do so within the context 
of the equilibrium tidal theory of Hut (1981), which is 
widely employed in the study of high-e migration (e.g., 
Eggleton et al. 1998; Wu& Murray 2003; Fabrycky & 
Tremaine 2007; Naoz et al. 2011). 

In §2 we summarize results from equilibrium tide the- 
ory, in the limit of weak dissipation, that are relevant 
for this work. We revisit the concept of the tidal quality 
factor Q and examine its relation to the time lag r of the 
equilibrium tidal response. The required r for high— e 
migration is estimated in §3. There we show that the re- 
quired strength of tidal dissipation for high-e migration 
is at least 10 times stronger than the inferred upper- 
limit for dissipation arising from the Jupiter-Io interac- 
tion. In §4 we briefly review some previous calculations 
of high— e extrasolar planet migration and comment on 
their prescriptions for tidal dissipation. We also discuss 
§5 some specific dynamical tide models in the literature 
and show that they cannot account for circularization of 
the hot Jupiters and thus are unlikely account for their 
presumed high-e migration. We summarize in §6. 

2. THEORY OF EQUILIBRIUM TIDES 

2.1. Physical Basis for Equilibrium Tides 

Hut (1981) asserts that the most attractive feature of 
the constant lag time model for tidal dissipation is its 
simplicity. That is, it is straightforward to implement in 
analyses of secular orbital dynamics. However, Socrates 
& Katz (2012) show that a constant single time lag r 
follows directly from the equations of motion, given the 
following basic assumptions 

• the tidal forcing and dissipation mechanism are 
slow and non-resonant 1 

1 Here, 'non-resonant' refers to the fluid response of the forced 
body. 
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• the equilibrium structure of the forced body is 
spherically symmetric 

• the dissipation is weak such that the decay time 
of the forced response is long in comparison to the 
characteristic forcing periods 

• only the quadrapolar component of the tidal force 
is important 

• non-linear effects are unimportant 

Altogether, these conditions lead to a constant lag time 
t of the forced response, which is an intrinsic property of 
the tidally forced object in question. Therefore, if Jupiter 
were to be placed in another system where the character- 
istic forcing frequency and amplitude are different than 
that felt by Jupiter from Io, the lag time r nevertheless 
remains the same as long as the assumptions above are 
correct. 

Socrates & Katz (2012) show that, given the assump- 
tions listed above, the constant lag time prescription fol- 
lows directly from the equations of motion for a spheri- 
cally symmetric forced fluid body. In the Appendix we 
show that, within the equilibrium tide approximation, 
quantification of the level of tidal dissipation with a tidal 
quality factor Q is only possible under the assumption of 
a single forcing frequency and a single constant lag time. 

2.2. relevant results of equilibrium tides 

Under the simplifying assumption of spin-orbit align- 
ment, the amount of energy dissipated per orbit is given 

by 



dt dE/dt 



(1) 



orb 



= E F (J^j x [A-2Bx + Cx 2 ] 

and the change in specific orbital angular momentum per 
orbit is 



AJ- 



orb 



= -j 



dt dJ/dt 
GM (M + )/a F , J 



(2) 



where E-p = 
^G(M + m pcr )a F , P F = 2n/y/G (M + m per ) /a| = 
2tt/tif and ap = a (l — e 2 ). Here x = fi/rtp where ft is 
the rotation rate and A — C are functions of e and given 
by eqs. A33 - A35 of Hut (1981). In the event of pseudo 
synchronous rotation where VL — VL VS — npB/C, orbital 
angular momentum J is conserved and ap is a constant 
and equal to the final circularization radius. Here M 
is the mass of the forced body, m per is the mass of the 
perturber and 

D ~ 3fc L rGm2 cr i?5 W 

is a characteristic time scale for tidal dissipation, where 
r is the lag time, k L is the Love number and R is the 
radius of the forced body. 

The evolutionary equations for the semi-major axis a 
and eccentricity e are obtained by dividing eqs. (1) and 



(2) by the orbital periods and then using Kepler's laws 
to relate the orbital energy E and J to e and a. The 
evolutionary equations in three limiting cases of interest 
are provided below. 

2.2.1. synchronization: e = 0, il ^ Q ps , M ^> m pcr 

Spin down of a major body by a minor body leads to 
expansion of the orbital separation between the two as 
a result of angular momentum conservation. The rate of 
this expansion is given by 



2(n-Q) M 



nt r 



(4) 



where n is the mean motion. The expression above de- 
termines the rate of orbital expansion e.g., of the Moon 
from the Earth. The ultimate source of energy and an- 
gular momentum originate from the spin of the major 
body. 

2.2.2. circularization: e -C 1, f2 = O ps , M <C m peT 

When e< 1, the pseudosynchonrous spin rate f2 ps w n 
and the rate of decay of orbital eccentricity follows 

i"t 

When the orbital angular momentum is large compared 
to the spin angular momentum, the forced body attains a 
pseudo synchronous state before circularization can pro- 
ceed. The source of energy for dissipation during migra- 
tion is primarily orbital. 

2.2.3. high-e migration: e — >• 1, ft = f2 ps; M -C m pcr 

As in the case of circularization the spin is pseudo syn- 
chronous during high-e migration and therefore, the or- 
bital angular momentum J remains a constant during 
inward migration. Consequently, the shape of the or- 
bit near periastron is approximately constant as long as 
e w 1 and therefore, the energy dissipated per orbit AE 
(e.g., Socrates et al. 2012) is approximately a constant 
as well. 

In the limit that e — > 1 the migration rate is given by 



a 
a 




(6) 



which is consistent with the expectation that AE w a 
constant per orbit. 

2.3. Q & t 

Since Goldreich & Soter's (1966) measurement of the 
tidal quality factor Q of Jupiter, Q has become a common 
parameterization for the level of tidal dissipation in gas 
giant planets and stars. In what follows, we discuss the 
relationship and applicability of Q and and the lag time 

T. 

In physics, the quality factor Q is often utilized to pa- 
rameterize the level of dissipation in a forced damped 
harmonic oscillator that obeys 



X 



to 2 X 



7 A = F(t) 



(7) 



.^oe - " + c.c. The fundamental frequency 



with F = ^\F^e lGt 

io and damping rate 7 are intrinsic properties of the os- 
cillator. When the forcing is slow such that uCw and 
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if the dissipation is weak, then the leading order solution 
is given by 



F(t) 



(8) 



with respect to the tidal problem, the expression above 
is analogous to the equilibrium tide. The term oc X is 
small in comparison to the X^ and does not contribute 
to energy transfer. To next order in forcing frequency 
a = 2n/P the contribution to X is due dissipation and 
is given by 



which allows us to write 

X= — — e lcrt - iS 

where the time lag is given by 



tF (9) 



+ c.c. (10) 



T = J/UJ 2 (11) 

and the phase lag 8 = or. 
For a damped, driven harmonic oscillator Q is defined 

as 



SdtdE/dt - AEc . 



2nE 



2ttE 



(12) 



where Ai? cyc is the energy dissipated per forcing cycle 
and is given by 



-AE,. 



dtXF ~ 2tt- 



(13) 



and the peak energy E , is the energy stored in the in- 
teraction 



E Q = (XF) = 



1 



dtXF = 



\F n 



and therefore, 



f dtXF 
J cyc 

' 2n(XF) 



8 = 



(14) 



(15) 



The oscillator Q is used as a proxy for the phase lag of 
the forced response. Its definition requires that the time- 
dependance of the forcing be described by a single Fourier 
harmonic. Given the forcing frequency a, a constant lag 
time r is also required to define Q. 

In general, the forced damped harmonic oscillator is 
not driven by a a single Fourier harmonic and the time- 
dependance can be completely arbitrary. The more gen- 
eral solution in the non-resonant weak friction limit may 
be written as 

X(t) = X^ +x « =jr<°> -rX(°) 

~l( 0) (t-r) (16) 

in the limit that r is small and thus the response X(i) 
follows the forcing F(t) by some time lag r. Note that 
the lag time r parameterizes the level of dissipation for 
a damped harmonic oscillator, for arbitrary forcing. The 
oscillator Q, which is equivalent to the phase lag 8, char- 
acterizes the level of dissipation for the specialized case 
of sinusoidal forcing. Unlike u and r, Q is not an intrin- 
sic property of the oscillator as it depends on a particular 



property of the external driving force i.e., the forcing fre- 
quency a. 

As previously mentioned, Socrates & Katz (2012) show 
that for the tidal problem, a constant lag time r is an 
intrinsic property of the forced body, given the physical 
requirements listed in §2.1, assuming the that fluid re- 
sponse of the forced body obeys the following equation 
of motion 

£+C-£ + D-£ = -VU T (17) 

where £ is the Lagrangian displacement and Ut is the 
tidal potential. The operators C and D can be thought 
of as a generalized "spring" coefficient and damping rate, 
respectively. 

The tidal potential can be decomposed into multipoles 
of spherical harmonic degree I. The lag time Tt of the 
corresponding tidal multipole response qi m is given by 



n = 



(18) 



where H^ m (x) — Vr £ Y£ m (fi) describes the spatial varia- 
tion of the tidal field about the co-ordinates, x = (r, ft), 
of the forced body. Both C and D are time-independent 
differential operators that depend, like the fluid density 
p, on the internal structure of the forced body. It then 
follows that, as in the case of the non-resonant weakly 
damped forced harmonic oscillator, the lag time is an in- 
trinsic property of a non-resonant weakly damped tidally 
forced body. 

Socrates & Katz (2012) also show that there is a con- 
stant time lag ti for each spherical harmonic degree I 
of the forced response. By adopting the approximation 
that the perturber is distant such that only the 1 = 2 
terms are kept, they are then able to reproduce the basic 
assumptions of Hut (1981). 

Even in the event that the time dependence of the tidal 
forcing is described by a single Fourier harmonic, it is 
not clear why, and under what conditions, the level of 
dissipation of the tidal problem can be parameterized by 
a quality factor Q. As previously mentioned, the tidal 
potential is the sum of the tidal harmonics, each of which 
may be a sum of Fourier harmonics. Furthermore, and 
on a more basic level, the left hand side of eq. 17 is not, 
in general, equivalent to the left hand of eq. 7. In other 
words, the tidal problem is not equivalent, in general, to 
a forced damped harmonic oscillator. 

Nevertheless, in the limiting cases where the constant 
time lag model is valid and when there is a single forcing 
frequency, a, the tidal quality factor Q is given by the ra- 
tio of the energy dissipated during the tidal forcing cycle 
and the average tidal interaction energy of the oscillating 
components give by (see Appendix) 

= (TT. (19) 



Q Eosc : fc2,m 2lT (llm^im) 

here qe m is the induced multipolar moment arising from 
the presence of the (normalized) tidal potential $ e m and 
the Eosc m denominator enforces the inclusion of only 
the oscillating components of qi m and ^em- Note the 
correspondence between the expression above and the 
expression for the oscillator Q given by eq. 15. In the 
Appendix we show that, under several limiting condi- 
tions, the tidal Q as defined above is equivalent to a 
single phase lag. 
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2.3.1. e=0; synchronization 

Both the m = ±2 components of the tidal potential os- 
cillate sinusoidally at the same frequency a = 2 |0 — n\. 
The to = perturbation is time-independent in this or- 
bital configuration and therefore, does not contribute to 
secular evolution. It follows that the tidal Q is well de- 
fined as discussed above. As we state in the Appendix, 
the peak energy Eq is due to equal contributions of both 
components. Furthermore, since the lag time r is inde- 
pendent of m, the forced response corresponding to each 
component of the tidal forcing possess the same phase 
lag (and angle) as well and therefore, the same Q. 

By using AE from eq. 1 and E from eq. All we may 
write the familiar relation 

— — = 2 \n-n\r. (20) 

In the expression above we made use of the fact the or- 
bital period and the period of the forcing cycle differ by 
a ratio of 2 \ fl — n\ jn in order to relate AE cyc to AE of 
eq. 1. 

2.3.2. e <C 1; circularization 

Again, the tidal Q is well defined for this orbital config- 
uration in part because the tidal forcing is described by a 
single forcing frequency a — n. The forcing is composed 
of three distinct Fourier harmonics that originate from 
the to = and to = ±2 components of the tidal poten- 
tial. As we state in the Appendix, the peak energy Eq is 
the sum of all three components, where the contributions 
from the to = and to = ±2 components are not equal. 
Furthermore, since the lag time r is independent of m, 
the forced response corresponding to each component of 
the tidal forcing possesses the same phase lag (and an- 
gle) as well and therefore, the same Q. Consequently, 
the level of dissipation can be parameterized by a single 
Q as well. 

By using AE from eq. 1 in the limit that = O ps (cf. 
eq. A 36 of Hut (1981)) and along with the expression 
for the peak energy E from eq. A14 we may write 

J-=nr. (21) 

circ 

Here, AE of eq. 1 is equal to AE cyc since the forcing 
frequency is equal to the mean motion n. 

2.3.3. e w 1; high-e migration 

In this limit, Q is not well defined i.e., there is neither 
a well-defined forcing cycle in order to compute an <f cyc - 
type integral nor a well-defined phase lag since the forcing 
is not sinusoidal. Nevertheless, it is possible to define a 
forcing cycle that is approximately periodic in the limit 
that e — > 1 since tidal potential and thus E vanishes at 
apoastron. In order to define a Q for high-e migration 
a peak energy E is needed, which we set to the peak 
interaction given by 

E = kL °™ R5 (22) 



p 



where r p = a(l 



e) is the periastron distance. Now we 
may define a tidal Q during high— e migration 

1 AE 3861 

7lF T. 



TABLE 1 

Inferred values of various tidal dissipation 
parameters for jljpiter-like planets a 





Jupiter-Io 


hot Jupiter 


high— e 


parameters 


interaction 


circularization 


migration" 


r (b) 

Q 


0.062 


0.25 


0.66 


5.9 x 10 4 


2.8 X 10 5 


1.1 x 10 5 


t D (Gyr) 


9.5 x 10 5 


20 


7.5 


T (hrs) 


1.2 x 10 4 


3.0 x 10 3 


1.1 x 10 3 


t v (hrs) 


1.3 x 10 5 


3.3 x 10 4 


1.2 x 10 4 


<7-l/ c (lCr 59 ) 


6.3 


25 


67 



a We take the Love number k L - 0.38, R p - 6.99 X lCTcm, 
M p - 1.9 X 10 30 g, m Io = 8.93 X 10 25 g, the mass of the stellar 
companion — Mq and the final circularization distance a p — 
0.06 AU. 

For high — e migration the tidal Q is not well-defined and we 
use Qhem as defined by eq. 23. 
c The parameter av has dimensions of (g cm 2 s) 1 . 



Interestingly, Q HEM as defined above remains approxi- 
mately a constant during high-e migration. Take for ex- 
ample, a case where the initial value of eccentricity is 
e = 0.999 and tidal dissipation shrinks the orbit to a 
value of e = 0.99. Though the orbital period changes 
by a factor of <~ 30, the value of Q HEM changes by less 
than 1%. The shape of the orbit near periastron - where 
all of the energy transfer takes place - remains approx- 
imately constant during evolution (e.g., Socrates ct al. 
2012). It follows that AE cyc = AE per orbit remains 
approximately a constant as does the peak energy E 
and therefore, so should Q HEM . 

2.4. other parameterizations 

Within the literature, there are other parameteriza- 
tions of the constant lag time for various theoretically- 
based motivations. Expressions relating them to the lag 
time t are given below. 

Hut (1981) defines a characteristic time scale for dissi- 
pation 

(24) 



T = 



GMt 



Eggleton et al. (1998; see also Fabrycky & Tremaine 
2007) use the viscous time ty 



t v =3(l + l/fc L )T = 3(l + l/fc L ) 



R 3 



(25) 



GMt 

and the coefficient ay incorporated by Hansen (2010; 
2012) 



°v = -k L rGR- 5 . 
3. APPLICATION TO GAS GIANT PLANETS 



(26) 



Q r 



2nE { 



o 



4096 



(23) 



In what follows, we determine the values of r, tv, and 
ay for Jupiter, based on its interaction with Io. We then 
compare these values with the required values for cir- 
cularization and high— e migration of Jupiter analogues 
(M = Mj, R = Rj and same structure). In particular, 
we estimate the value of Q and to for these cases. 
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3.1. Jupiter-Io interaction 

Jupiter's rotation is super-synchronous with respect 
to its orbit with Io (fij > n la ). In order to obtain a 
timescale for the evolution of Io's orbit, use eq. 4 and 
integrate from a = to a = a\ Q in order to write 



a . 13 

~ h \a=a lo = 2A*0 

where t Q — 4.5 Gyr is the age of the Sun and 



Ganymede 



A = 



i— Io 



Lio 



4.3 



(27) 



(28) 



where Li is the orbital angular momentum if the i th 
Galilean moon. In obtaining an upper limit for dissipa- 
tion, the inclusion of the factor A conservatively assumes 
that the current resonant configuration of the inner three 
Galilean moons is preserved during outward migration. 
From the form of to given by eq. 3, we obtain a limit of 
the lag time r J Io for the Jupiter-Io interaction 

Tj-jo < 0.062 s (29) 

where we have set = 0.38. By using eq. 20 we may 
determine the lower-limit for Q,_ Io , the tidal Q for the 
Jupiter-Io interaction 

~>4 



<9j_ Io > 5.9 x 10 4 



(30) 



as found by other authors (Goldreich & Soter 1966; Yo- 
der & Peale 1981; Leconte et al. 2010). 

3.2. circularization of hot Jupiters 

Nearly all hot Jupiters with orbital periods P < 5 days 
and semi-major axes a < 0.06 AU have circular orbits, 
while the values of e belong to a broad distribution for 
gas giants at larger separation. A common interpreta- 
tion is that planets that are relatively close-in circularize 
quickly due to tidal dissipation's strong dependance on 
orbital separation. Consequently, it is possible to obtain 
a separate constraint on the dissipation parameters r, Q, 
etc., by presuming that tidal dissipation is responsible for 
circularizing Jupiter-like planets with P < 5 days within 
a characteristic age to = 10 Gyrs of these systems. 

If we set t^ 1 = e/e, given by the right hand side of eq. 
5 then we can write an expression for to 

t D =2t . (31) 

With the help of eq. 3 and by setting fc^ = 0.38 and 
ap = 0.06 AU we arrive at a constraint for r 

t> 0.25 s (32) 

for a Jupiter-analogue in orbit around a Sun-like star. In 
the above constraint, we chose = 0.38 and ap = 0.06 
AU. Given this value of r, eq. 21 enables us to determine 

Q 

Q circ < 2.8 x 10 5 . 



(33) 



3.3. high-e migration of hot Jupiters 

A leading theory for the formation of hot Jupiters 
is high-e migration (see Socrates et al. 2012 for a re- 
view). In particular Fabrycky & Tremaine (2007) pre- 
dicted that, as a result of high-e migration, the orbital 
angular momentum of hot Jupiters should, in general, 
be mis-aligned with the spin axis of the stellar host, 



which was later confirmed observationally by Winn et 
al. (2010). 

The migration time t m between an initial eccentricity 
ei and some final value eo can be found by numerically 
integrating eq. 6 as depicted in Figure 1 of Socrates et 
al. (2011). For example, a Jupiter analogue starting at 
an orbital period P ~ 12 yrs that circularizes to a final 
distance of ~ 0.06 AU with eo ~ 0.1 migrates in a time 
t m given by 

t m = t(e A = 0.994) - t(e = 0.1) « ^t D . (34) 

For a fixed value of J or of, the migration time t m does 
not vary considerably if the initial value of the orbital 
period varies at the order unity level. For high— e mi- 
gration to operate the above expression requires that to 
satisfies 



or the lag time obeys 



tD<^to 



t > 0.66 s 



(35) 



(36) 



for a Jupiter analogue. 

With our definition given by eq. 23 along with the 
value of r above to place a limit on Q during high-e 
migration 

Q HEM < 1.1 x 10 5 . (37) 

Note that t m and consequently, r and Q HEM are rather 
insensitive to a. 

3.4. discussion 

If indeed Jupiter analogues with P F — 5 days success- 
fully undergo high— e migration, then according to the 
results of §§3.1 and 3.3, there is at least an order of mag- 
nitude disagreement between the two requirements. In 
other words, the time r required for high— e migration 
is ten times larger than the allowed upper-limit for r of 
Jupiter. For example, if Jupiter's lag time r J _ Io = 0.66 
s, as required for high-e migration, then its tidal qual- 
ity factor would decrease by an order of magnitude such 
that Q,_ Io < 6 x 10 3 . 

The tension between two requirements for r is obscured 
when tidal dissipation is quantified by Q. At face value, 
the requirements for Q from circularization and high— e 
migration indicate that less energy is dissipated per cycle 
than that inferred from the Jupiter-Io calibration. How- 
ever, as detailed in §2.3, the tidal Q, unlike the lag time 
r, is not an intrinsic property of the forced body in the 
equilibrium theory. Consequently, if the determination 
of t from the Jupiter-Io interaction is correct, then Ta- 
ble 1 reveals that a Jupiter analogue would not be able 
to migrate from its initial position to Pp m 5 days within 
the Hubble time. 

Note that the requirement on r from tidal circular- 
ization of hot Jupiters is marginally consistent with the 
value inferred from the Jupiter-Io interaction, which is 
roughly in agreement with several previous studies (e.g., 
Leconte et al. 2010; Hansen 2010 & 2012). 

4. PREVIOUS WORK ON HIGH-ECCENTRICITY 
MIGRATION 

The fact that the lag time r necessary to migrate a 
Jupitier analogue to a final circularization period Pf ~ 5 
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days is at least 10 times larger than tj_i does not 
present an enormous hurdle for high-e migration scenar- 
ios. Due to the strong dependence of tp on the final 
circularization radius of and planet radius R p , it is pos- 
sible that a factor of 10 may be recovered in order to 
make the migration time i m consistent with to < 10 Gyr. 

However, there is a general understanding in the liter- 
ature that for modest values of tidal dissipation inferred 
from the Jupitcr-Io interaction (<5 t _ Io ~ a few x 10 5 ), 
the migration time of a Jupiter analogue with Pp w 5 
days is of order - 300 Myrs (Wu & Murray 2003; Fab- 
rycky & Tremaine 2007) in stark contrast with our results 
of the last section. 

4.1. Wu & Murray 2003 

Wu & Murray (2003; hereafter WM) produced the first 
calculations of high-e migration for hot Jupiters, with 
particular emphasis on HD 80606b. 2 Their evolutionary 
equations are taken from Egglcton et al. (1998). In this 
model, tidal dissipation is ultimately modeled after the 
orbit-averaged constant time lag model of Hut (1981). In 
order to determine the dissipation rate, a constant time 
lag t must be specified. 

However, inspection of their eq. (A9) reveals that the 
value of r in their calculations is 



1 



T WM 



192. 



3 x 10 5 
Qo 



P 



12yrs 



s (38) 



where n is the mean motion and they set Qo = 3 x 10 . 
Note that r WM is not a constant during evolution and 
may be ~ 3000 times larger than upper limit for t t _ Io 
given by Table 1. This particular deformation of the 
constant r model of equilibrium tides results in an infinite 
amount of energy dissipated per orbit in the limit that 
e — > 1 and P — > oo for a fixed periastron. Finally, in this 
e -> 1 limit, the tidal Q HEM = -2irE /AE 0. 

Similar errors can be found in related works (Wu et al. 
2007; Wu & Lithwick 2010). For example, consider eq. 
5 of Wu et al. (2007), which they state is equivalent to 
eq. 10 of Hut (1981), in fact reveals that 



d\nt 



dt 



dint 



Wu 



dt 



P 
x — . 

Hut Pf 



(39) 



Since the orbit is pseudo synchronized during high— e 
migration, Pp, J ad of are all fixed during evolution. 
Therefore, for a orbit starting with P = 12 yrs migrat- 
ing on a track of angular momentum J corresponding to 
Pf = 5 days, the lhs of the expression above exceeds the 
rhs by a factor ~ P/Pf ~ 10 3 . 

4.2. Fabrycky & Tremaine 2007 

Like WM, Fabrycky & Tremaine (2007; hereafter FT) 
employ the evolutionary equations of Eggleton et al. 
(1998), where the prescription for tidal dissipation is 
based upon Hut (1981). Unlike WM however, they did 
not alter the orbit evolutionary equations. They qualita- 
tively reproduce the presumed migration history of HD 

2 In a planetary context, high— e migration was first utilized to 
explain the misalignment of Triton's orbit with the spin of Neptune 
by Goldreich et al. (1989). Those authors do not provide justifica- 
tion for their prescription for tidal dissipation. In the high— e limit, 
the behavior of their evolutionary equations is, correctly, identical 
to that of Hut (1981). 



80606b calculated by WM. In order to do so, they used 
a viscous time ty of 

t VFT = 0.001 yr = 8.8 hrs. (40) 

The qualitative agreement between WM and FT seems 
to strengthen the case for high-e migration with respect 
to the history of HD 80606b. When computing the dis- 
tribution of the inclined hot Jupiters, FT use a value for 

ty of 

ty FT = 0.01 yr = 88 hrs. (41) 

However, while FT used the correct equations of mo- 
tion their value for ty departs from the value determined 
from the Jupiter-Io interaction given in Table 1 

t Vj _ IO = 1.3 x 10 5 hrs, (42) 

by four orders of magnitude. Furthermore, note that the 
value FT use for ty in both cases leads to a value of 
\AE/E \ ~ 0.1 during high-e migration. 

In addition, FT write an expression in their equation 
A10 for the tidal Q that follows 

Q oc P (43) 

which as previously discussed at length, is incorrect 
within the context of high-e migration. 

4.3. discussion of previous work 

The mis-calibrations highlighted above demonstrate 
the danger of using Q when parameterizing the strength 
of tidal dissipation. For the case of the synchronization 
tide in the Jupiter-Io system, there is little quantitative 
difference between using Q and r since the forcing pe- 
riod of Jupiter is approximately constant through the 
evolution of Io's orbit. Consequently, in the equations of 
motion, replacing lag time - or phase lag - in the equa- 
tions of motion can be done without loss of accuracy e.g., 
as in Goldreich & Soter (1966). Nevertheless, due to the 
manner in which Q is defined, its placement within equa- 
tion of motion (cf. §11 of Goldreich & Soter 1966) can be 
extremely confusing or incorrect altogether. 

The intuition derived from Goldreich & Soter (1966) 
for the synchronization tide of Jupiter cannot be applied 
to the more general tidal problem, where the orbit is 
eccentric. In the case of the Jupiter-Io interaction, a 
phase lag and physical lag angle are equal to one another, 
both of which are equivalent to the energy dissipated over 
a well-defined periodic sinusoidal forcing cycle. There 
is, in general, no well-defined forcing cycle in the tidal 
problem and therefore, no single phase lag. Furthermore, 
for the example of a pseudo synchronized eccentric orbit, 
the lag angle is a function of orbital phase, where its 
value at e.g., apoastron and periastron are not equal. 

Despite previous mis-calibrations in the level of tidal 
dissipation in migrating super-eccentric gas giants, as 
previously mentioned, Table 1 indicates that a modest 
increase by a factor of <~ 10 in the relative strength of 
dissipation allows for Jupiter analogues to migrate from 
a period of 12 years to 5 days within the age of the 
Universe. This is a direct consequence of the fact that 
tidal dissipation is an extremely strong function of or- 
bital angular momentum J, periastron r p or a F . For ex- 
ample, a F = 0.071 AU in Figure 1 of FT, rather than HD 
80606b's value for a F = 0.061 AU, dissipation is also in- 
trinsically weaker since they chose M p = 7.8 Mj, migra- 
tion was completed in ~ 3 — 4 Gyrs, rather than to = 10 
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Gyrs, at long orbital periods the fractional amount of 
time migrating was small since the planet underwent 
Kozai oscillations and therefore, spent a significant frac- 
tion of its history at high— J. By accounting for all these 
effects, well over two of the orders that separate eqs. 40 
and 42 may be erased . In what follows, we discuss some 
of the consequences for studies of gas giants that result 
from the requirement of modestly enhanced dissipation 
for high-e migration. 

5. SOME IMPLICATIONS FOR STUDIES OF GAS GIANTS 

5.1. tidal theory 

It was realized early on that turbulent convection is 
severely inadequate in providing the required amount of 
dissipation in order to account for Jupiter's tidal Q (Gol- 
dreich & Nicholson 1977). Since then, the leading can- 
didate for the source of tidal dissipation has been the 
resonant excitation of low frequency, nearly incompress- 
ible, waves that are primarily restored by the Coriolis 
force. This topic is technically challenging and continues 
to be in a state of development (Ogilvie & Lin 2004; Ar- 
ras 2004; Wu 2005a,b; Goodman & Lackner 2009; Ogilvie 
2009). 

The results and relations of equilibrium tide theory 
cannot be directly compared with predictions of dynam- 
ical tidal theories. Nevetheless, the value of AE re- 
quired for circularization for a given Pp within the Hub- 
ble time is largely independent of any particular tidal 
model. Goodman & Lackner (2009) write an expression 
for a tidal Q = —AE/2ttEo during hot Jupiter circu- 
larization that results from the inertial wave model. In 
their expression for Q, their peak energy Eq is the peak 
energy in the fluctuating component of the interaction 
energy, appropriate for circularization given by our eq. 
A14 in the Appendix. 

Their reported value for Q is Q > 3 x 10 7 (e.g., Good- 
man & Lackner 2009) for Jupiter analogues (core radius 
R c < 0.2Rj) with Pp = 5 days. Table 1 indicates that 
the required Q for circularization of hot Jupiters with 
Pp w 5 days is Q < 3 x 10 5 . Thus, calculations of tidal 
circularization of gas giant planets arising from the in- 
ertial wave mechanism fall short by a factor of 100 or 
so. High-e migration will then be almost certainly more 
difficult to accomplish by the inertial wave mechanism. 
However, a calculation of the forced response that leads 
to the inertial wave mechanism for a highly eccentric or- 
bit has yet to be performed, so a quantitative comparison 
with results from equilibrium tidal theory are therefore, 
not yet possible. 

A clear difference between the properties of the tidal 
potential felt by gas giants undergoing high-e migration 
and Jupiter is the amplitude of the tidal potential. The 
ratio is given by U la /U± ~ 10~ 3 . If indeed super-eccentric 
migrating Jupiters experience enhanced tidal dissipation 
during migration, then it may be due to non-linear ef- 
fects. 

5.2. the inflated radii problem 

The hot Jupiters are inflated in that, unlike Jupiter, 
their radii are, in general, significantly larger than the 
zero temperature solution. There are several solutions 
to this problem, many of which depend upon tidal dissi- 
pation in one way or the other. 



There is a class of mechanisms that can be thought 
of "delayed" contraction mechanisms. Here, the planet 
is brought into its close orbit in such a way that its de- 
generacy is lifted (Burrows et al. 2007; Leconte et al. 
2010; Ibgui & Burrows 2009; Batygin & Stevenson 2010; 
see also Wu & Lithwick 2011). One candidate is tidal 
dissipation due to circularization during high-e migra- 
tion. Table 1 indicates that in order for these mecha- 
nisms to work, tidal dissipation must be much stronger 
for super-eccentric migrating Jupiters in comparison to 
that inferred from the Jupiter-Io interaction. 

Thermal tidal torques can render the hot Jupiters in a 
persistent state of asynchronous spin. The resulting on- 
going dissipation of the gravitational tide, presumably at 
great depth, allows for a steady source of thermal energy 
required to inflate the radii significantly above the zero 
temperature solution (Arras & Socrates 2009a, 2009b; 
2010). In order for this mechanism to work, the out-of- 
phase quadrupole induced by thermal forcing near the 
optical photosphere must be equal to the out-of-phase 
quadrupole induced by dissipation of the gravitational 
tide. Such an equilibrium can only be realized if the 
relative level of tidal dissipation in the hot Jupiters is 
comparable to that inferred from the Jupiter-Io interac- 
tion. Therefore, if high-e migration is indeed responsible 
for delivering gas giants from initial orbital periods of 
w 12 yrs to a final circularized state of « 5 days, then 
thermal tidal torques cannot compete against the gravi- 
tational tidal torques. All of this, however, is under the 
assumption that the lag time for the synchronization and 
circularization tide are equal, which may not be the case 
in reality. 

6. SUMMARY 

If high-e migration leads to the production of hot 
Jupiters with final circularized orbital periods Pp w 5 
days, tidal dissipation in these extrasolar gas giants must 
be at least 10 times stronger than that inferred from the 
Jupiter-Io interaction. While this alone cannot rule out 
high— e migration, it is in severe conflict with several pre- 
vious studies of high-e migration. We showed that the 
discrepancy between our analysis and these studies arises 
from the fact that they employed dissipation strengths 
which are 10 3 — 10 4 larger than the inferred upper limit 
of the Jupiter-Io interaction. 

While current theories of tidal excitation and dissipa- 
tion in gas giant planets can account for the inferred level 
of dissipation in the Jupiter-Io interaction (Ogilvie & Lin 
2004) they cannot currently account for the amount of 
dissipation required for the observed circularization of 
hot Jupiters. If high— e migration is in fact responsible 
for hot Jupiter migration, the discrepancy between the 
level of computed dissipation in the inertial wave mech- 
anism and that required to account for migration is even 
larger. If, for example, thermal emission from migrat- 
ing Jupiters in orbit about late-type stellar hosts are ob- 
served by on-going and future direct-imaging campaigns 
(Dong et al. 2012), then the tidal dissipation time to 
cannot be much larger than to J$ Gyr. 

In the Appendix we point out that the the tidal quality 
factor is well-defined and can be related to a phase lag in 
only in two limiting cases: asynchronous and circular or- 
bits or pseudo-synchronous orbits with small eccentricity. 
We find the relation between the tidal quality factor Q 
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and the lag time r for these two limiting cases previously 
mentioned and additionally during high-e migration. 

Given the results presented here, there is now tension 
between various sub-fields of gas giant studies: the theory 
of tidal friction in gas giant planets, few-body dynamics 
of exoplanet systems, theory of inflated radii, the evolu- 
tionary scenario outlined by Goldrcich & Soter (1966) for 
the Jupitcr-Io interaction and even perhaps planet for- 
mation. Either the standard thinking in one (or more) of 
these categories requires revision or there are extra ingre- 
dients at work in forming (stellar spin-orbit) mis-aligned 
hot Jupiters. The production of hot Jupiters is not a 
solved problem, though high-e migration is still viable. 
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APPENDIX 



INTERACTION POTENTIAL, LOVE NUMBER, PEAK ENERGY AND ENERGY TRANSFER RATE 

We review the mathematical machinery necessary to address the tidal problem. We assume the forced object is 
perturbed by a source of gravity of mass m at distance D(t) in a Keplerian orbit with true anomaly $(£) and that the 
strength of the tidal field is weak. Therefore, only terms that are responsible for linear perturbing forces and linear 
responses are kept. The main purpose of what follows is to provide the necessary physical ingredients that are required 
to connect the results of equilibrium tidal theories that utilize a constant time lag model for tidal dissipation (Hut 
1981) with those models that utilize a tidal quality factor Q (cf. Goldreich & Soter 1966). 

interaction potential 

To leading order in perturbation theory, the interaction potential is given by (cf. Newcomb 1962) 

Ri = J (PxpZ-VUt (Al) 

where J d 3 x is taken over the volume of the forced body, p(x) is its static fluid density and £(x, t) is the Lagrangian 
displacement field. The tidal potential Ut is given by 

U T = -Gra £ ^L. Y e * m (tt/2, $) Y tm («) 

= -S^( n )*<mW ( A2 ) 

lm 
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which serves to define ^t m {f). Note that we we have specialized to the case where the perturber is located at co- latitude 
9 = 7r/2. Now, we may write the interaction energy as 

Hj = -Y,4m^tm (A3) 

where qe m is the multipole moment (cf. Press & Teukolsky 1977) 

q em (t) = I d 3 xp £ • Vr e Y* n («) . (A4) 
Note that for the £ = 2, the ^ m 's carry dimensions of frequency squared. 

ke - the Love number 

The Love number ke is defined as 

UUR) = h U<t(R) (A5) 

where U* nd is the additional gravitational potential resulting for the induced tidal response of degree I and R is the 
un-perturbed radius of the forced body. The induced potential is given by (e.g., Jackson 1999) 

4-7T 

C/ind = -GV — — qtmr-^Ytm (fi) (A6) 

em M + L 

and with this we can solve for each qe m in terms of ke, *£ m and R 

qe m = ^em/we ( A7 ) 

where 

2 _ 47rG 

Wt = fc £ (2^+l)^ +1 - ( 8) 

£o - the peak energy 

The tidal Q is a dimensionless number that parameterizes the relative level of dissipation of the forced body over 
some interval in time. Given a tidal theory, a solution of the equations of motion allows for a determination of the 
energy dissipated AE over a forcing cycle. In order to construct the dimensionless parameter Q, which is equivalent 
to a phase lag, a characteristic energy Eq is needed. Typically, E is referred to as the "peak energy." It is equal to 
the average energy in the tidal interaction or the peak energy of the potential energy of the flow. When computing 
Eq only fluctuating components of the tidal potential are considered in order to maintain the standard relationship 
between Q and the phase lag. In the literature, Eq is defined in many ways and for the sake of clarity, we provide a 
definition for it. 

Insert the expression for the induced multipoles in terms of k and * em from eq. A7 into the expression for Hi and 
obtain 

H i = -l^ k '—^ q = -2J**"l l w t- ( A9 ) 

em em 

The interaction energy Hi has contributions from both time-independent and fluctuating components. We define the 
peak energy Eq with the following expression 

E ^-(Hi) = ^I dt ]T |*L|M 2 = E l*L|M 2 (Aio) 

eye J eye OSCj £ m oscjm 

where the integral is taken over the forcing cycle and only terms resulting from the oscillating component of the tidal 
potential are kept. 

synchronization: e = 0, Q ^ fl ps 

Tidal synchronization for circular orbits is due to forcing from the I = 2 and m — ±2 portions of the tidal potential, 
which are responsible for its time- varying component. In this case, the oscillating component of the tidal field can be 
thought of as a being due to a single Fourier harmonic with forcing frequency a = 2 \ n — fi|. The peak energy in the 
interaction Eq is then given by 

E = -(Hi)= ]T \Vem\ 2 /™ 2 e 

e=2,m=±2 

_3k e Gm 2 R 5 

Eo ~- tf— (AU) 

During synchronization, there is an m = tidal deformation as well. However, since it is docs not oscillate, it cannot 
contribute to energy transfer and consequently, does not factor into the definition of Eq. 
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circularization: e<l, !! = f2 ps 

For slightly eccentric orbits with e < 1 in pseudo synchronous spin states with Own, isolation of the oscillating 
terms of Ut can be accomplished by a Fourier decomposition. In the rotating frame, let 



( 



a\ eim <t, = j2 Xi m {e)e l(yk ^ t (A12) 



DJ 

k= — oo 



where the forcing frequency dkm ~ kn — mft and X^J n (e) are the Hansen coefficients, which are given by 
Xl m {e) = ^fj' n dte^-^-t 

S km + I [{t + 1 + 2m) 6 k , m+1 + ((. + ! -2m) <5 fe , m -i] + O (e 2 ) (A13) 

to leading order in eccentricity. In the e -C 1, CI = Cl ps w n limit, the only terms that contribute to time-dependent 
driving are, for I = 2, the (m = 0, fc = 1), (m = 2, fc = 2 ± 1) and the (m = —2, fc = — 2 ± 1) terms. Note that all of 
the oscillation frequencies are therefore given by a^ m = ±n. 
With decomposition into Hansen coefficients, the interaction energy Hj and peak energy Eo becomes 

£ = -<#!>= l^trnkf/lXe 
osc, £mk 

^npqr^ (Al4) 

where the ^tmk's are defined by the relation 

= GmW em ^ = Gm Z tm x i m (ey akmt = E *™ ^ (A15) 
u k a k 

and the Fourier transform of the induced multipole moments qt m k are 

<Um ^Y.^^ 1 ( A16 ) 

k 

each of which satisfy 

qimk = ^emk/t&e- (A17) 

high-e migration: e — > 1, = ri ps 

For the case of high-e migration, the time dependance of the tidal potential is far from that of a single Fourier 
harmonic. Its magnitude is highly peaked at periastron and falls of quickly with increasing distance. We choose to 
define the peak energy as the entire interaction potential, since all of the Yi m 's vary with time and therefore, contribute 
to the potential energy available to perform work on the orbit. We write 

E Q = max [-Hi] = k 2 (A18) 
where r p = a(l — e) is the periastron distance. 

energy transfer rate 

The rate of energy transfer E between the orbit and internal degrees of freedom of the forced body is given by 

E = - J (PxpZ-VUt- (A19) 
The definition of Ut in terms of the ^e m 's may be employed in order to write 

E = -52qt m 9 tm . (A20) 

CONDITIONS FOR A SINGLE Q FROM MANY OSCILLATORS AND THE TIDAL PROBLEM 

Even for the case of semi-diurnal forcing (I — 2; m = ±2) during synchronization of a body in a circular orbit, the 
forcing may be thought of as originating from two driving terms. It is not clear that the level of tidal dissipation in 
such a system can be summarized by a single Q, which itself is composed of a single phase lag. In order to proceed, a 
definition for Q is required for a set of decoupled forced damped harmonic oscillators, each of which obey 

X n + ulX n + ln X n = F n (t) = ^ e* 7 "* + c.c. (Bl) 
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where w„ is the fundamental frequency of each oscillator amplitude X n , j n is the corresponding damping rate and a„ 
is the driving frequency. Each of the n oscillators has a Q = Q n given by 



1 _ l yc dE n /dt -i yc X n F n 



S n (B2) 



Qn 2nE ,n 27r(X„F„) 

where r„ = 7n/w 2 is the lag time and S n is the phase lag. In the expression above X n (t) is the solution in the 
non-resonant weak friction limit where X n (t) — F n (t — t„)/oj^. 
Now, define Q for the entire set of oscillators 

1 ErJ^dEJdt = E n j yc dE n /dt = E„MF , K | 2 /^ 

Q En^Eo.n E„27T(X„F„) £JF ,„| 2 /^ 

The expression above is well-defined only if there is a single forcing cycle and therefore, only a single forcing frequency. 
In the event that there is a single lag time for each oscillator, then 5 n — a constant. It follows that 

^ = Sn = S. (B4) 

Or in other words, a single Q that is equivalent to a single phase lag quantifies the level of dissipation of a set decoupled 
oscillators only if their lag times are equal to one another. 
To make contact with the tidal problem, we restrict to I — 2 and write the tidal Q as 



1_ iyJE/dt _§ cyc dE/dt _ f cyc dE/dt £*= 2 , ro £ yc <7| ro ^ 



Q 2nE 2tt (if,) £ £=2>ro 2tt (<G m * tm ) E OS c,fc2, m ^ tf m *<m ' 



(B5) 



As in the case of a set of decoupled forced non-resonant weakly damped harmonic oscillators, a single Q for the 
combined system only makes sense if there is a single tidal forcing frequency. The time derivative q* lm enforces that 
the Ylim m t ne numerator runs over only the oscillating components of the tidal potential and the corresponding tidal 
response. Expand into Fourier harmonics and write 



j_ _ J2e=2,rnk S Jrnk IjVmfcj 

Q £osc,^=2, mk I ^imk I l w "l 

Each individual phase lag 5e m k — <JmkTt are equal to one another in the event that their respective lag times as equal 
to one another as well. If this is the case, we may then write 

q = 5f.rnk = S. (B7) 

Comparison of eq. B6 with eq. B3 shows a correspondence between a set of non-resonantly forced, weakly damped 
harmonic oscillators (X ni F n ) with the tidal problem (qe m ,^£m) in the equilibrium weak-friction limit. Though the 
general tidal problem is not equivalent to a set of decoupled harmonic oscillators with constant time lags, their 
respective expressions for Q are equivalent under these limited conditions. That is, a single tidal Q that represents a 
single phase lag only holds when e<l and ft = n for circularization or e = during synchronization. In addition, 
the time lag for each of the m components of the tidal force of degree I = 2 must be equal to one another. From this, 
we arrive at the following conclusion: a single Q can serve to parameterize the level of tidal dissipation only if there 
is a single forcing frequency a and a single constant lag time r. That is, parameterization of tidal dissipation with a 
constant r is applicable in a more general class of orbital conditions in comparison to parameterization with a tidal Q. 



